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1 Introduction a tuned comparison at partonic and hadronic levels 



The current theoretical status of Drell-Yan [1] physics 
is widely overviewed in the resent papers [2,3] where 
the necessity of further in-depth study of the Drell- 
Yan-like (DY) processes is emphasized. 

At hadron-hadron colliders DY processes serve for 
the normalization purposes [4, 5] , provide information 
about weak interactions [6,7], contribute to the back- 
ground in many searches for new physics beyond the 
Standard Model (SM). At the LHC they will be used 
also for precision fits of partonic density functions in 
the regions of x and factorization scale values which This article is the second ste P in the series of SANC 



at one-loop precision was realized between four codes: 
[12], [15], [16] and [23,24]. 

To the contrary, a detailed tuned comparison for 
neutral current DY NLO either at parton or at hadron 
levels is still on the way. In Ref. [25] one can find a 
comparison between the results of the HORACE and 
ZGRADE codes for the inclusive NC DY cross section 
at the hadronic level. Partial results of the comparison 
at the partonic level between SANC and another code 
can be found in [26] . 



were not yet accessed experimentally. P a P ers devoted to DY processes. First, we presented 

The theoretical calculations of the DY processes in Rcf - I 27 ] thc CC casc - Here we show results for the 

for high energy hadronic colliders were performed at NC case - Moreover, in thc DY branches of SANC we 

the level of one-loop QED and electroweak (EW) ra- can stud y the interplay between QCD and EW NLO 

diative corrections (RC) by several groups, see papers corrections within thc same framework. Thc current 

[8-17] and references therein. QCD corrections are status of this stud y was presented in the report to 

known up to thc next-to-next-to-leading order [18-20]. ATLAS MC working group at CERN [28]. Detailed 

For thc experiment, both charge current (CC) and description of thc QCD branch for DY CC and NC 

neutral current (NC) DY processes represent great in- processes will be presensed elsewhere, though the first 

terest. Both types of the processes can be easily de- results arc alre ady published in [29]. 

tected. CC processes have a larger cross section, but The paper is organized as follows. In the second 

NC processes have a clean dimuon signature and pro- section we demonstrate the implementation of the EW 

vide an additional information, e.g. about the weak NLO corrections into the SANC framework [23]. The 

mixing angle. location of the standard set of the SANC modules: 

The theoretical study for CC DY is pushed ahead FF (form factors), HA (helicity amplitudes), BR and 

more then for NC — several independent codes exist MC (brcmsstrahlung), is shown on the screenshot of 

for EW NLO RC, see [3]. Within thc 2005 Lcs Houchcs the SANC Processes in the EW 4f NC sector. Ibi- 

workshop [21] and the 2006 TEV4LHC workshop [22] dem we show the rather short analytical expression 
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for the differential cross section of the hard photon an event generator, based on the FOAM package [30] . 

contribution. The generator itself will be described elsewhere (see [31]) 

SANC team has the advantage of experience in and made accessible from the project homepages after 
the calculation of brcmsstrahlung. We can obtain the intensive rolling and testing, 
numerical results both by the semi-analytical expres- 
sions (it is the standard BR module of the SANC) 2 Neutral current Drell-Yan processes 

and by a Monte Carlo integrator or generator (with 

2.1 Born level 

help of the new MC module). 

At first we will consider interactions of free quarks 

In the third section we describe the adopted "sub- 
bartons). At the leading order (LO) the unpolarized 

traction" procedure in the expressions for virtual, soft 

cross section of the partonic subprocess q{pi)+q{p2) — * 

and hard contributions to the NC DY cross section. 

(7, Z) -> £(p 3 ) + £(p 4 ), q = (u, d, c, s, b) is given by 

Numerical results are presented in the Sect. 4 both 
for the Born and for the EWNLORC. We investigated 0{ - > 9§ P ' 



V (s) + ^V A (s) 



(1) 



the independence of the result on the variation of the where 

soft-hard separator u> which provides one of the most . , -,2 

s = -{Pi+P2), 



important checks of the calculation. The sensitivity of / 4m? 

(3 s = Jl-^-, (2) 



s 

the corrections to the variation of "subtracted" quark 

and rai being the lepton mass. Here we denoted 

mass singularities was also studied. 

In the conclusion we overview the present status of Vo,a(s) — Q q Qi + 2Q q Qiv q veRex(s) 

implementation of NC DY calculations into the SANC + + («/ ± 4) I x(») | 2 > 

system. Vq = 1^ - 2Q q sin 2 O w , 

As the main point of this article we offer for the Vg = _ 2Q i sin 2 9w , 

import a stand-alone code for NC DY EW NLO RC _ T ( 3 ) 

a q — 1 q 1 

at the partonic level together with the environment in (3) 

a e =iy. (3) 

which it was run. The sketchy description of this code 

The Z/7 propagator ratio Xz{s) with constant (or s- 

is presented in Appendix. Codes are accessible from 

dependent) Z width reads 

SANC project homepages http : //sane . j inr . ru and 

1 s 

http://pcphsanc.cern.ch. Xz{s) = - - - Tj^jr - ( 4 ) 

For production of numbers at the hadronic level where s w and c w are the sine and cosine of the weak 
the SANC team created a Monte Carlo integrator and mixing angle Qw 
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2.2 EW Radiative Corrections at the Partonic Level 

As the rule of the SANC basement, we subdivide the 
EW RC into the virtual (loop) ones, the ones due to 
soft photon emission, and the ones due to hard photon 
emission. Later on in the section 4 in the Tables 1 
and 2, we demonstrate the independence of the EW 
RC off an auxiliary parameter u> which subdivides the 
soft and hard photonic contributions. 

In Fig. 1 we show the location of the 2/ — * 2/ NC 
processes at the SANC tree. 



? 



Root 

\t SANC 

S QED 

t <a bn 

» □ Precomputation 
<9 Processes 
*■ □ 3 legs 
3 4 legs 

t ^3 Neutral Current 
o- □ Bhabha 

■a ti ti -> f f 

Fmf1 f1 -> f f (FF) 
Fmfl fl -> f f (HA) 

Fwf1 f1 ->f f (BR) 

rmfl fl -> f f (MC) 
Q- □ Charged Current 
*■ □ 2f2b 
□ QCD 



Fig. 1. SANC tree for the process 2/ -> 2/ 

Moving along the menu sequence SANC — > EW 
— » Processes — > 4 legs — * 4f — * Neutral Current 
->fl fl -> f f (FF,HA,BR,MC) user arrives at the 
first standard SANC module, the scalar Form Fac- 
tors (FF); then at the second module, Helicity Ampli- 
tudes (HA); then at the third module, the integrated 



bremsstrahlung (BR); and finally at the fourth mod- 
ule, the fully differential bremsstrahlung (MC). 

• FF and HA 

Earlier, in our article [23] we presented the Covari- 
ant (CA) and Helicity Amplitudes (HA) of fxfiff — > 
NC process, with all 4-momenta being incoming, for 
any of its cross channels s,t or u. The expressions for 
the CA and HA (see Eq.(30) and (33) of the last ref- 
erence) of this process can be written in terms of six 
FF. 

• BR and MC 

The BR module computes the soft and inclusive 
hard real photon emission: 



q(pi) + «(P2) -» l(pa) + Kpa) + 7(P5), 



(5) 



where the momenta of corresponding particles are given 
in brackets. We do not discuss the soft photon contri- 
bution here, referring the reader to the system itself. 
As far as hard photons are concerned, we realized two 
possibilities of the integration over its phase space: the 
semi-analytical one (BR) and the one by means of a 
Monte Carlo integrator or generator (MC). 

The first one is based on two different sets of ana- 
lytical integrals over two choices of kinematic variables 
parametrizing the hard photon phase space. In the 
first option (CalcScheme=0), the phase space looks 
like 



2tt 



where the two-body phase spaces are: 



(6) 



(1) 1 VA(a,y,0) l^ 

d<P\ 1 — dcosu 

1 8tt s 2 



7 ' 
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1 y/X(s',mj,mj) 1 



(2) 



-d cos 3 d(j) 3 . (7) 



8ir §' 47r 

Here, s' = — (p 3 + pi) 2 and 9 1 is the cms angle of the 
photon. Since, for the chosen parametrization of the 
phase-space the process matrix element squared does 
not depend on the angle </> 7 - cms azimuthal angle of 
the photon - it is already integrated out in the phase- 
space. The angles 9 3 and (j) 3 define the orientation of 
the momentum p 3 in the rest frame of the compound 
(P3 + P4=0). These parameters vary in the limits: 



with the Born cross section given by Eq. (1). The cou- 
pling functions are: 

Ai(s, s') = 2a q a e Re[ X (s') - x(s)] 
A 2 (s, s') = 2Q q Q g a q a t Re[ X {s') + *(*)] (11) 
+8v g a g veatRe [x(s)x* («')] ■ 

We give also a more simple expression neglecting lep- 
ton masses. Then all masses remains only in the argu- 
ments of logarithms: 

S2 1 s/2 





a 


ds' 


ir 


X 




'in J- 






. m 


1 a 3 






3 .s 2 





Q„Qi I Re 



< 0j,6 3 < it, 0<(f) 3 <2ir, 
\m\ <s'<§- 2Vsw . (8) 

After the integration over all angular variables, we ob- 
tained the compact expression for a single differential . , f\ * f'^ 
distribution of the hard photon contribution to the 
NC DY process, dohard/ds'. We neglected terms pro- 
portional to the initial quark mass but kept terms pro 
portional to the lepton mass: 



Q 2 l n 1 a (s') +Qj In— - 1 a (s) 



Q q Qe{xz(s) +Xz{s') 



(12) 



do-hard _ a s 2 + s' 2 
ds' ir s 2 s- 



Q 2 ln-^-1 Uo(s') 



m„ 



with cto(s) being the massless limit of Eq. (1). 

In the second option (CalcScheme=l), the module 
BR computes the distribution: d<Thard/dc, where c = 
cos(Zp2P4). The phase space in this case reads: 
1 



+Q 2 (L' - (3 S ,) -a (s) 



d$^ 



dcds 1 dZ±d§ 1 , 



(13) 



2a 3 

T 



H — —QqQe ~ 3 



£+ 



Ms,s') f^|±Lj 



-08 



+Q q QiA 1 (s 1 s')^L' 
4a 3 2 m 2 

x{v (s) 



{ 

-V A (S) 



, s_ s — m 
2 I — - 2 

s s_ 



where 



4 ^ + 3^ \L' B - 3^ft 



La = In 

I -0s 



s± = §± s' . 



2 9 7T 4 S 

where Z4 = — 2p 4 p 5 and 7 is the cms azimuthal angle 
of P5 varying within 27r limits. After integration over 
(j> 7 , Z4, varying within limits 



^_(l-ft)<Z 4 <^_(l+ft), (14) 



and s' — within the same limits as in Eq. (8), one 
gets very cumbersome result for the single differential 
distribution over variable c, where one has to neglect 
terms proportional to all masses keeping them only in 
(10) the arguments of logarithms. The user may see this re- 



(9) 
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suit after corresponding run of the SANC BR module 
with CalcSchemc=l. 

The MC module provides fully differential hard 
bremsstrahlung contribution to the partonic cross sec- 
tion. The contribution is given in a form suitable for 
further numerical integration or simulation of events 
in a Monte Carlo generator. 

3 Treatment of quark mass singularities 
3.1 Partonic level 

To perform the subtraction procedure at the partonic 
level cross section 



0"! — <5"0 + OSV + O'hard , 



(15) 



we proceed in the same way as in our previous paper 
on DY CC [27]. 

The subtracted expression Z\<r MS from the com- 
plete calculations with massive quarks: 



Aa 



MS 



-21n(l-&)-l 



i-e 

<7o(s&) 



where 



Ami 



(16) 



(17) 



Next, Qi and rrn are the charge and the mass of the 
given quark; M is the factorization scale; <7o(s£i) is 
the cross section at the partonic level with the reduced 
value of s — > s£j. 



The plus prescription in Eq. (16) can be treated in 
the following way: 



Aaf s = lim (Aa sv + Aa hald 



MS 



(18) 



The first contribution, Aasv, related to soft and 
virtual photonic contributions, is given by 



A*sv = -Qi In — » In 



M 2 , 4w 2 1, 2 4w 2 



■ In" 



s 2 s 



_l n _ + _l n _ + 2) (To ( S ), (19) 



where we took into account that mi = m 2 = m q and 
Qi = Q% = Q\ for the NC case. 

The second one, Z\<7h a rd, related to hard photon 
emission, is 



Aa] 



hard 



£mi„ 



1 + ^2 M 2 

1 — r ln — 



where 



-21n(l-0-l ao(*0. 



s — 2y/soj 



(20) 



(21) 



Using subtraction procedure, the cross section with 
O (a) corrections is given by 



•MS ' a ' J\ 

a 1 = o\ — Zla 1 



MS 



(22) 



Then it can be convoluted with PDF's in the usual 
way. An equivalent subtraction procedure applied at 
the hadronic level will be discussed right below. 

3.2 Hadronic level 

The differential cross section of the DY process at the 
hadronic level can be obtained from the convolution 
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of the partonic cross section with the quark density 
functions: 



1 1 



dc 



9192 



xq 2 (x 2 ,M 2 ) 1 ' 
dc 



any kind of an experimental cut. When the returning 
to the Z-resonance is allowed by kinematic cuts, the 
natural choice of the factorization scale is M 2 = M|. 
For the region of higher invariant masses of the lepton 
JO{c,xi,x 2 ), (23) pair, Mj+j-, it is better to take M 2 



J J dxi dx 2 qi(x u M 2 ) 



where the step function 0(c, x\,x 2 ) defines the phase 
space domain corresponding to the given event selec- 
tion procedure. The partonic cross section is taken 
in the center-of-mass reference frame of the initial 
quarks, where the cosine of the muon scattering an- 
gle c is defined. The transformation into the angle c 
defined in the cms of an initial hadrons involves the 
Jacobian: 



dc Ax x x 2 

J = tt = 5— > a = xi + x 2 - c(xi - x 2 ), 

dc a z 



l-(l-c) 



2xi 



SX\X 2 . 



(24) 



The parton densities with bars in Eq. (23) mean the 
ones modified by the subtraction of the quark mass 
singularities: 



1 — z \ ra z q 



(25) 



where q(x, M 2 ) can be taken directly from the exist- 
ing PDF's in the MS scheme (see Ref. [10] for the 
corresponding formula in the DIS scheme). It can be 
shown analytically, that this procedure is equivalent 
to the subtraction from the cross section, given by 
Eq. (22). In the approach with subtraction from PDF's 
it is easy to keep the completely differential form of 
the sub-process cross section and therefore to impose 



4 Numerical Results 

The input parameters set is taken the same as used in 
Ref. [21] (leshw_input .h, see Appendix): 



G F 


= 1.16637 x 10- 5 GeV~ 2 , 






a(0) 


= 1/137.03599911, 




= 0.1187, 


M w 


= 80.425 GeV, 


r w 


= 2.124 GeV, 


M z 


= 91.1867 GeV, 


r z 


= 2.4952 GeV, 


M H 


= 150 GeV, 


m t 


= 174.17 GeV, 


m e 


= 0.51099892- 10- 3 GeV, 


m u 


= TTtd — 66 MeV 




= 0.105658369 GeV, 


m c 


= 1.55 GeV, 


m s 


= 150 MeV, 






m T 


= 1.77699 GeV, 


m b 


= 4.5 GeV, 


\v ud \ 


= \V CS \ =0.975, 


\v us 


= \V cd \ = 0.222. 



(26) 



4.1 Partonic level 

We begin with a presentation of several numerical re- 
sults derived at the partonic level, investigating the 
independence of the results of some unphysical param- 
eters and on the choice of an EW scheme. First, we 
investigate the independence off the parameter CD, and 
the residual dependence on the initial quark masses 
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after the procedure of subtraction of quark mass sin- dence arises, in principle, which becomes stronger if 
gularities (see section 3) by choosing for quark masses cuts get tighter. The numbers in the Tables 1 and 2 
two values: first as in Eq. 26, and the second — ten are derived for the only cut M^+^- > 50 GeV. As seen, 
times lower. the residual quark mass dependence is very weak. 

All the numbers were obtained with the aid of 
4.1.1 Independence off parameter u> and quark masses a FORTRAN code which consists of a hand-written 



The sum of the "soft" and "hard" photon contribu- 
tions to the total and differential cross sections should 
not depend on the soft-hard separator u>, which wc 
varied from 1(T 3 ^ GeV to 1(T 5 ^ GeV. We ob- 
served that for u> = 1(T 4 ^ GeV and u> = 1CT 5 ^ 



main and FORTRAN modules automatically gener- 
ated by s2n . f software of the SANC system, see Ap- 
pendix. 

4.1.2 a(0), G F , G' F schemes 



GeV the numbers for the one-loop corrected cross sec- Here we study the EW scheme dependence of the cor- 
tions &i agree within four digits shown, and therefore rected cross section a 1 arising from the definition of 
we present only one row in the Tables 1 and 2. For EW constants in the a(0), G F and G F schemes, 
the relative radiative correction factors there is a tiny In Table 3 we show the results for a and <7i at the 
dependence on lu, and we present the RC for both w's Born and one-loop levels and for the corresponding 
using marks 1) for 1(T 4 and 2) for 1(T 5 . The differ- RC factors S in three EW schemes, 
ence between the results for two w's is, however, much The a(0) and G F EW schemes are defined as usu- 
below any reasonable estimate of the theoretical accu- ally [32] . In the G' F scheme one assigns the same one- 
racy, loop value of the coupling constant standing at all pho- 
The cross sections are shown in picobarn and the ton vertices a QE D ~ 1/132.544. It has been adopted 
radiative corrections in percent. All the results in Ta- in Ref. [12] and used in Rcf. [27] for the sake of com- 
blcs 1 and 2 are derived in a(0) EW scheme and after parison only. 

the subtraction of quark mass singularities in the WS The results for the radiative corrections are derived 

scheme. in MS subtraction scheme, and the factorization scale 

Each Table contains four rows of 5, two pairs for is taken to beM z . There are notable deviations of the 

different ui for m u = m d = 0.066 GeV and m u = m d = corrected cross sections from the corresponding Born 

0.0066 GeV. After subtraction of quark mass singular- values and between the Born values themselves in the 

ities the total cross section should be fully independent three schemes. One can observe a certain degree of 

off quark masses. If cuts arc applied, a residual depen- stabilization of one-loop corrected cross sections as 
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Table 1. The total lowest-order and one-loop corrected cross sections fro and fri for the process uu — > /i + /i~(7) 
in the a EW scheme and corresponding relative one-loop correction S for Cu = 10~ 4 -t^ GeV, rows 1) and 
oj = lCT 5 ^! GeV, rows 2). Cut value: M M+M - > 50 GeV. 



V§, GeV 


70 


90 


110 


500 


1000 


2000 


5000 


14000 


ao, pb 


4.312 


369.6 


4.995 


7.324xl(T 2 


1.806xl(T 2 


4.499xl(T 3 


7.191xl(T 4 


9.171xl(T 5 


fri, pb 


4.713 


397.2 


4.732 


8.403X1CT 2 


2.105X1CT 2 


5.256xlCT 3 


8.269xlCT 4 


l.OlOxlCT 4 


m q = m u 


s, %, 1) 


9.303 


7.468 


-5.259 


14.74 


16.61 


16.84 


14.99 


10.09 


5, %, 2) 


9.303 


7.481 


-5.261 


14.74 


16.61 


16.84 


14.99 


10.09 


m q = m u /10 


<5, %, 1) 


9.296 


7.466 


-5.261 


14.74 


16.61 


16.84 


14.99 


10.09 


«5, %, 2) 


9.296 


7.480 


-5.263 


14.74 


16.61 


16.84 


14.99 


10.09 



Table 2. The total lowest-order and one-loop corrected cross sections ao and &i for the process dd — » (7) 
in the a EW scheme and corresponding relative one-loop correction 8 for lj = 10~ 4 -^ GeV, rows 1) and 
Co = lQ- 5 ^- GeV, rows 2). Cut value: M fl+fl - > 50 GeV. 



V§, GeV 


70 


90 


110 


500 


1000 


2000 


5000 


14000 


fro, pb 


2.895 


472.9 


5.237 


3.968xl0~ 2 


9.601xl0~ 3 


2.381xl0~ 3 


3.802xl0~ 4 


4.847xl0~ 5 


fri, pb 


3.079 


497.3 


5.318 


4.439xl0~ 2 


1.090xlO~ 2 


2.692xl0~ 3 


4.149xl0~ 4 


4.863xl0~ 5 


m q = m d 


8, %, 1) 


6.328 


5.164 


1.551 


11.88 


13.53 


13.04 


9.142 


0.319 


S, %, 2) 


6.327 


5.167 


1.550 


11.88 


13.53 


13.04 


9.141 


0.317 


m q = md/10 


5, %, 1) 


6.329 


5.162 


1.552 


11.88 


13.53 


13.04 


9.142 


0.318 


5, %, 2) 


6.329 


5.165 


1.551 


11.88 


13.53 


13.04 


9.141 


0.317 
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Table 3. The total lowest-order parton cross section ao and one-loop corrected cross section &\ in pb and relative 
RC factor 8 for the two EW schemes Gf and G' F for u) — 10~ 5 -^ GeV and invariant mass cut M M + M - = 50 GeV; 
for details see the text. 



Vs, GeV 


70 


90 


110 


500 


1000 


2000 


5000 


14000 


uu^> n + n (7) 


Q 


A ,% 


-0.55 


-1.41 


-0.34 


-0.59 


-0.66 


-0.64 


-0.53 


-0.30 


Gf 


a , pb 


4.434 


394.9 


5.259 


7.464xl0~ 2 


1.839xl0~ 2 


4.581xl0~ 3 


7.323xl0~ 4 


9.338xl0~ 5 


ai, pb 


4.739 


402.9 


4.748 


8.453xl0~ 2 


2.119xl0~ 2 


5.290xlO~ 3 


8.313xl0~ 4 


1.013xl0~ 4 


8, % 


6.878 


2.019 


-9.718 


13.25 


15.23 


15.47 


13.53 


8.443 


G' F 


ao, pb 


4.608 


395.0 


5.338 


7.826xl0~ 2 


1.929xl0~ 2 


4.807xlO~ 3 


7.685xl0~ 4 


9.800xl0~ 5 


ai, pb 


4.805 


403.2 


4.759 


8.657xl0~ 2 


2.176xl0~ 2 


5.442xl0~ 3 


8.576xl0~ 3 


1.048xlO~ 4 


8, % 


4.285 


2.075 


-10.83 


10.61 


12.76 


13.21 


11.60 


6.913 


A, % 


+ 1.01 


+0.07 


+0.23 


+2.41 


+2.69 


+2.87 


+3.16 


+3.46 


dd — > (7) 


a 




-0.96 


-1.29 


-1.02 


-1.06 


-1.09 


-1.07 


-0.91 


-0.045 


Gf 


ao, pb 


3.055 


505.3 


5.573 


4.144xl0~ 2 


1.002xlO~ 2 


2.485xl0~ 3 


3.966xl0~ 4 


5.057xl0" 5 


di, pb 


3.109 


503.8 


5.373 


4.487xl0~ 2 


1.102xlO~ 2 


2.721xl0~ 3 


4.187xl0~ 4 


4.885xl0~ 5 


8, % 


1.781 


-0.306 


-3.593 


8.268 


9.994 


9.509 


5.550 


-3.408 


G' F 


ao, pb 


3.094 


505.3 


5.596 


4.240xlO~ 2 


1.026xlO~ 2 


2.545xl0~ 3 


4.062xlO~ 4 


5.18OX10" 5 


<ti, pb 


3.126 


503.7 


5.371 


4.534xl0~ 2 


1.115xl0~ 2 


2.756xl0~ 3 


4.240xlO~ 4 


4.940xl0~ 5 


8, % 


1.026 


-0.310 


-4.027 


6.932 


8.718 


8.291 


4.368 


-4.619 


A, % 


+0.55 


-0.02 


-0.04 


+1.05 


+1.18 


+1.29 


1.27 


+ 1.13 
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can be seen from the rows of A values, showing devi- scheme. Besides input parameters Eq. (26), one has 
ations of o\ in a(0) and G' F schemes from <7i in Gf to specify cuts on the lcpton invariant mass M M + M -, 
scheme: muons transverse momenta pt(m) and muons pseudo- 

rapidity r](fi) which were used at production of figures 



A(a, G'p) = *!(a, G' F )/<j\{G F ) - 1,% . (27) 

However, for the process uu — > /U + pT (7) , espe- 
cially at high energies, the G' F scheme does not seem 
to work satisfactorily. 



below: 



Af„ + „- > 50GeV, 
Pt (a* ± ) > 25GcV, 



Note, that for a(0) scheme it is sufficient to show |f?GO| < 1.2. (28) 

only Z\-rows, since all the other entries are already 



given in Tables 1 and 2. 

As concerning the choice between the a(0) and Gf 
schemes, we favor the latter since its effective energy 
scale is known to be of the order of a weak boson mass. 
Moreover, in the results presented in the a(0) scheme, 
the hadronic contribution to the vacuum polarization 
in the photon propagator is taken into account in an 
approximation, using the one-loop formulae supplied 
by the quark masses listed in Eq. (26). In any case 
the difference between the schemes is related to the 
contribution of higher orders of the perturbation the- 
ory and gives us a crude estimate of the corresponding 
contribution to the theoretical uncertainty. 

4.2 Hadronic level 



In Figs. 2 and 3 we show the distributions da/dpr 
and dS/dpr obtained with the aid of the MC integra- 
tor, while in the Figs. 4 and 5 the same distributions 
obtained with the aid of the generator. The corre- 
sponding distributions agree with each other within 
statistical errors, which are larger for the case of the 
MC event generator. 

In the Figures 6 and 7 we show the distributions 
da/dM and d5/dM over invariant mass M — M^+^- 
for the case of the MC integrator, and in the Figures 8 
and 9 — for the case of the generator. 

One may draw the same conclusions as for the dis- 
tributions in px- 

The similar distributions are presented in Ref. [17] 
calculated, however, with different setup and cuts, the- 
refore they cannot be compared straightforwardly. 



5 Conclusions 



Results presented in this section were obtained with 
the aid of two standalone packages: 1) a MC integrator 
based on Vegas algorithm [33] and 2) a MC generator 

based on FOAM algorithm [30]. All the calculations In addition to the results presented here, one has to 
are done with Les Houches 2005 setup [21] inGj? EW take into account one more O (a) contribution of the 
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Distribution of u* transverse momentum 




$5 30 35 40 45 50 



Fig. 2. Muon transverse momentum distribution for yfs = 
14 TeV, as obtained from the MC integrator. Both the 
Born (black line) and the one-loop results (red line) are 
shown. 



Distribution of u* transverse momentum 




WeV] 



Fig. 3. Relative corrections 8 as a function of the muon 
transverse momentum pr(At + ), as obtained from the MC 
integrator. 



Distribution of u* transverse momentum 



> 




$5 30 35 40 45 50 



P T r)[GeV] 

Fig. 4. Muon transverse momentum distribution for y^i = 
14 TeV, as obtained from the MC event generator. Both 
the Born (red line) and the one-loop results (blue line) are 
shown. 



Distribution of u* transverse momentum 




25 30 35 40 45 50 



WeV] 

Fig. 5. Relative corrections 5 as a function of the muon 
transverse momentum pr{^ + ), as obtained from the MC 
event generator. 
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Distribution of invariant mass of u*u' pair 



Distribution of invariant mass of u*u' pair 




100 110 120 130 140 150 
M(nV) [GeV] 




%0 60 70 



100 110 120 130 140 150 
M(nV) [GeV] 



Fig. 6. Invariant mass distribution of the M + A t_ pair at Fig. 8. Invariant mass distribution of the A l+ A t_ pair at 
y/s — 14 TeV, as obtained from the MC integrator. Both yfs = 14 TeV, as obtained from the MC event generator, 
the Born (black line) and the one-loop results (red line) Both the Born (red line) and the one-loop results (blue 



are shown. 



line) are shown. 



Distribution of invariant mass of iiV pair 




50 60 70 



90 100 110 120 130 140 150 
M(nV) [GeV] 



Distribution of invariant mass of u*u' pair 



£ 100 

10 



-20 






Lii 



50 60 70 80 90 100 110 120 130 140 150 

M(nV) [GeV] 



Fig. 7. Relative corrections 5 as a function of the invariant Fig. 9. Relative corrections S as a function of the invariant 
mass M M + M -, as obtained from the MC integrator. 



mass M„+„- , as obtained from the MC event generator. 



photon-induced subprocess 7 + q — * q + I + I, which 
was recently computed in Ref . [34] . 



In this way the complete set of one-loop electroweak 
radiative corrections to neutral current Drell-Yan pro- 
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cesses has been computed within the SANC system u, 14 — d, 16 — fi~, etc. The package sanc_dy_nc_vl . 00 
environment. Our results are implemented into semi- is accessible for the following set of particles: 1313_1212, 
analytical FORTRAN codes as well as into a Monte Carlo 1313.1616, 1414.1212, 1414.1616. 
event generator. The codes can help to increase an Total set of hies inside the package is: 
accuracy of the theoretical description of the SM pro- Instruction hies : 
cesses, which is required for the forthcoming LHC data 
analysis. Tuned comparison of our results for NC DY 
with the ones of analogous calculations performed by 
other groups, started for the CC case in Refs. [21,22], 
is in progress within the Les Houches 2007 workshop. 
The comparison should help us to derive a common 
conclusion on the resulting theoretical uncertainty in 
our predictions. s2n_declare . h 



README 

RELEASE-NOTES 
CHANGES 
LICENSE . TXT 
INSTALL 

Declaration hies: 



Initialization and various input hies: 



Appendix 



s2n_init . f 
sanc_input .h 
leshw_input . fa- 
te v41hcw_ input .h 



Libraries for various functions , including Vegas inte- 
gration, see INSTALL hie in the package. 
Main hies: main_nc_vegas _xx_yy . F 
SSFM originated from 



In the Appendix we present a Technical Description 
of the sanc_dy_nc_vl . 00 package intended for calcu- 
lation of the total DY NC cross section at the partonic 
level. 

The main aim of this description is to demonstrate 

how SANC Standard FORTRAN Modules (SSFM) 

can be used by other codes. We also give a short guide 

to the main flags which are being used not only in this 

package but also govern the work of any package which nc_f f _xx_yy . F (FF) 

uses SSFM, as our integrators and generators do. nc_si_xx_yy . f (HA) 

• Overview of the package nc_br_xx_yy . f (BR) 

All hies are produced by s2n package of SANC [this file contains three SSFM (subroutines) 

project (v. 1.10) except main*_xx_yy.F, declaration hies nc_bo_xx_yy (...), nc_br_xx_yy (...), 

— *.h and libraries — *.a. Here "xx" and "yy" stand nc_ha_xx_xx_lspr (...)] 

for the standard SANC held indices: 12 — e~, 13 — nc_ha_xx_yy . f (MC) 
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The steps of calculation in the main* files in accor- 
dance with Eqs. (22), (18), (15) are: 

- step of declaration and initialization 

- step born is realized by flag iborn=l, 

(To is computed by integration over c of Eq. (13) of 
the function_lc_xx_yy, 

— via SSFM call nc_br_xx_yy (...,born,...) 

- step asv is realized by flag iborn=0, 

it is computed by integration over c of 
the function_lc_xx_yy, 

— virt via SSFM call nc_si_xx_yy (...,sigma) , 
(inside this module there exists a call to SSFM 
nc_ff_xx_yy (...)) 

— soft via SSFM call nc_br_xx_yy (...,soft,...) 

- step Aa^ 

by calculation of Ajgy through Eq.(19) 

- ^ep Aa^ ld 

by integration over £ of the function_ls_xx_yy 

— via SSFM call nc_bo_xx_yy (...,bornk), 
see Eq. (20) 

- step do-hard/dg' 

by integration over s' of the function_lspr_xx_yy 

— via SSFM call nc_ha_Kx_yy_lspr (...,hard) 
or alternatively 

- step rfCThard/^ (3) 

by integration over 4d-phase space of Eq. (6)-(7) 
of the function_4d_xx_yy 

— via SSFM call nc_ha_xx_yy (...,hard) 

• Options of the flags 
iqed(I) - choice of calculations for QED correction: 
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1=0 without QED correction 
1=1 with all QED correction 
1=2 with ISR QED correction 
1=3 with IFI QED correction 
1=4 with FSR QED correction 

iew(I) choice of calculations for EW correction: 
1=0 without EW correction 
1=1 with EW correction 

iborn(I) - choice of scheme of calculations of the par- 
tonic cross section: 

1=1 only Born level 

1=0 Born + 1-loop virtual corrections 

gfscheme(I) - choice of the EW scheme: 
1=0 a(0)scheme 
1=1 Gf scheme 
1=2 G F scheme 

ilin(I) - choice of the linearization at the calculation 
of the partonic cross section: 

1=0 without linearization 
1=1 with linearization, i.e. neglecting 
spurious terms 0{a 2 ) 

ifgg(I) - choice of calculations of photonic vacuum 
polarization T gg : 

I=-l — 

1= — 1 

1= 1 — 1 + LFgg 

1= 2 - 1/(1 - kF gg ) 

with k = - 9 ,> . 

16V 
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ihard(I) - types of the hard bremsstrahhmg phase- 
space integrations: 

1=1 integration over s' 
1=4 4d integration 

isetup(I) - choice of the setup: 
1=0 Standard SANC 
1=1 Les Houches Workshop 
1=2 TcV4LHC Workshop 

The package can be accessed from project home- 
pages http:/ /pcphsanc. cern. ch/ download/ sanc_ dy_nc 
_vl.00.tgz or http://sanc.jinr.ru/download/sanc_dy 
_nc_vl.OO.tgz. 
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